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ABSTRACT 



We report the discovery of a giant, loop-like stellar structure around the edge- 
on spiral galaxy NGC 4013. This arcing feature extends 6' (~ 26 kpc in projected 
distance) northeast from the center and 3' (~ 12 kpc) from the disk plane; likely 
related features are also apparent on the southwest side of the disk, extending 
to 4' (~ 17 kpc). The detection of this low surface-brightness (/zr = 27.0lo;2 
mag/ ) structure is independently confirmed in three separate datasets from 
three different telescopes. 

Although its true three dimensional geometry is unknown, the sky- projected 
morphology of this structure displays a match with the theoretical predictions for 
the edge-on, projected view of a stellar tidal streams of a dwarf satellite moving 
in a low inclined (~ 25°), nearly circular orbit. Using the recent model of the 
Monoceros tidal stream in the Milky Way by Pefiarrubia et al. as template, we 
find that the progenitor system may have been a galaxy with an initial mass 
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6 X 10^ Mq, of which current position and final fate is unknown. According to 
this simulation, the tidal stream may be approximately ~ 2.8 Gyr of age. 

Our results demonstrate that NGC4013, previously considered a prototypical 
isolated disk galaxy in spite of having one of the most prominent Hi warps de- 
tected thus far, may have in fact suffered a recent minor merger. This discovery 
highlights that undisturbed disks at high surface brightness levels in the optical 
but warped in H I maps may in fact reveal complex signatures of recent accretion 
events in deep photometric surveys. 

Subject headings: galaxies: individual (NGC 4013) — galaxies: dwarf — galaxies: 
evolution — galaxies: interactions — galaxies: disk galaxies — galaxies: warps 

1. INTRODUCTION 

In the last decade, the study of the formation and evolution of the Milky Way (MW) 
has been revolutionized by the first generation of wide-field, digital imaging surveys. The re- 
sulting extensive photometric databases have revealed for the first time the existence of spec- 
tacular stellar tidal streams (e.g., that from the Sagittarius dwarf galaxy; Ibata ct al. 2001; 
Martfnez-Delgado et al. 2001; Majewski et al. 2003) as well as large stellar sub-structures 
in the halo (Ncwberg et al. 2002; Rocha-Pinto ct al. 2004; Juric et al. 2008), interpreted as 
the fossils of the hierarchical formation of our Galaxy. The discovery of the Monoceros tidal 
stream (Newberg et al. 2002; Yanny et al. 2003), located close to the Galactic plane outside 
the MW disk (as well as similar structures seen around the M31 disk, Ibata et al. 2005), 
indicates that mergers might play a relevant role in the formation of the outer regions of 
spiral disks (Peiiarrubia et al. 2006) like that of the MW. Moreover, inside-out disk forma- 
tion from continual accretion of in-falling material is an observed feature in cold dark matter 
(ACDM) simulations of the growth of structures on galaxy- sized scales (e.g., Abadi et al. 
2003), and is now also a common feature of galactic chemical evolution models (e.g., Alibes, 
Labay, & Canal 2001; Chiappini, Matteucci, & Romano 2001) attempting to explain trends 
in disk chemical abundance patterns. These various results provide clear evidence that the 
destruction of satellite galaxies plays a relevant role not only in the formation of MW-hke 
spiral galaxies generally, but for their disks as well as their halos. Furthermore, these results 
suggest that the stellar mass assembly of the MW disk, and disks in general, likely continues 
actively to the present epoch. 

The accretion of satellite galaxies on low-inclination orbits may be an important for- 
mation mechanism of galactic disks in the ACDM paradigm (e.g., Abadi et al. 2003). The 
existence of low-inchned tidal streams in the outer regions of our Galaxy (e.g., the Mono- 
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ceros or Triangulum/ Andromeda tidal streams, Yanny et al. 2003; Majewski et al. 2003; 
Grillmair 2006) supports this prediction, although the external origin of these over-densities 
is still a matter of debate, mainly due to the severe extinction hindering the exploration 
of low-latitude areas. Recently, for example, Kazantzidis et al. (2008) and Younger et al. 
(2008) have suggested that, in a ACDM context, the bombardment by dark matter sub- 
halos of an initially cold disk formed at z ~ 1 would also yield the formation of low-surface 
brightness substructures {/i ^ 25 mag/ □" ) with filamentary shape that may locally resemble 
tidal streams in phase-space. Clearly, additional observational data (e.g. detailed chemical 
abundances) obtained along the Monoceros stream are needed to distinguish between both 
theoretical scenarios. 

The search for extragalactic analogies to MW minor mergers is required not only to (1) 
show that the MW is not unusual in this respect, but to (2) obtain externally-viewed "snap- 
shots" of different phases of such interactions, (3) explore the range of possible mass/orbit 
combinations for such activity, and (4) estimate the fractional contribution of accreted mass 
and the mass spectrum of such events in the life of MW-like systems, an issue that remains 
unresolved (Majewski 1999). In this way, a systematic survey of tidal streams around other, 
nearby disk galaxies can provide new constraints and insights on the hierarchical formation 
and structure of MW-like galaxies beyond the previously limited views afforded by our own 
Galaxy and M31 (Ibata et al. 2007). A need to build up statistical information on the number 
and distribution of tidal streams in galaxies is driven by the availability of state-of-the-art, 
high resolution cosmological simulations that offer, for the the first time, the opportunity 
to use these observations to probe the theoretical predictions of galaxy formation in the 
framework of the ACDM paradigm (e.g.. Bullock & Johnston 2005). 

Promising galaxies in the hunt for extragalactic tidal streams are those displaying out- 
standing asymmetries in optical or Hi images. It has been long suggested that these per- 
turbations are a result of gravitational interaction with nearby companions. In some cases, 
apparently isolated galaxies exhibit morphologies more commonly associated with interact- 
ing systems. Recently, the prototypical isolated, warped disk galaxy, NGC 5907, was found 
to be surrounded by a spectacular stellar tidal stream (Shang et al. 1998; Martfnez-Delgado 
et al. 2008), but one with no obvious companion (which may have been completely dis- 
rupted). In the case of both the MW and M31 there are also obvious disk warps — both 
stellar and gaseous — but no immediately obvious perturber, though small, nearby, likely 
tidally disrupting companion satellites have been suggested as a cause in each case — (Sato 
& Sawa 1986; Ibata & Razoumov 1998; Baihn 2003; Bailin & Steinmetz 2004). These ex- 
amples suggest that spirals with warped disks but no large nearby companions, may have 
undergone minor mergers in the last Gyrs. 
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Under this prepense, we have been observing disk galaxies to faint surface brightnesses in 
the search for evidence of minor mergers. We (Martmez-Delgado et al. 2008) have previously 
reported our observations of the debris of a minor merger around the NGC 5907 disk system. 
Here we report similar work on the edge-on disk, NGC 4013. NGC 4013 is one of the 62 
luminous {Mb < —16.9) members of the Ursa Major cluster of galaxies, a nearby, late-type 
dominated, and low mass galaxy cluster. According to the HYPERLEDA database, the mean 
hehocentric radial velocity of NGC 4013 (836 kms~^) places it at a distance of 14.6 Mpc. 
It is also a relatively isolated system, with the two nearest, slightly more luminous, cluster 
neighbours (NGC 4051 and NGC 3938) at ~ 170/250 kpc projected distance away. Classified 
as an Sb galaxy with a maximum observed rotational velocity of 195 kms~^ (Bottema 1996), 
an extinction corrected total absolute magnitude of -20.1 Mabs sband (Vcrhcijen & Sancisi 
2001), an optical scale-length of 2.8 kpc (40") (van der Kruit & Searle 1982) and an isophote 
25 mag arcsec~^ radius of 5.2'(de Vaucouleurs et al. 1991) , NGC 4013 is very similar to 
the Milky Way. NGC 4013 is, moreover, famous for its prodigiously warped Hi disk, with a 
line of nodes close to parallel with that of the line of sight and with one of the largest warp 
angles observed (~ 25deg). On one side the warp extends out to about 8 kpc (2') from the 
nominal plane (Bottema et al. 1987; Bottema 1995, 1996). The box/peanut shaped bulge 
of NGC 4013 with its 'X'-hke morphology indicates the presence of a bar being observed 
edge-on (Patsis & Xilouris 2006). 

Here we report the optical detection of a faint, low-surface brightness, loop-hke struc- 
ture that appears to be part of a giant, low- inclination stellar tidal stream of a disrupted 
dwarf sateUite looking very similar to the expected form of the Monoceros tidal stream. 
This looping structure suggests the likelihood of there having been a previous interaction of 
NGC 4013 with a low-mass companion. Given the fact that NGC4013's very prominently 
warped H I layer matches in shape and orientation those of the apparent tidal loop makes this 
galaxy system a very compelling example of the likely link between disk warps and mergers, 
and an interesting case study for models of warp formation. 



2. OBSERVATIONS AND DATA REDUCTION 



We observed NGC 4013 with three different telescopes and instruments. In each case 
the same structures appeared, albeit with different degrees of clarity. 
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2.1. KPNO 0.9-m Telescope 

NGC4013 was initially observed with the Kitt Peak National Observatory (KPNO) 
(now WIYN) 0.9-m, f/8 telescope as part of a pilot survey of low surface brightness features 
(like stellar tidal streams) around warped, nearby disk galaxies by CP. and S.R.M. The 
initial sample included NGC3044, NGC3079, NGC3432 and NGC4013, all edge-on systems 
of large angular size with Hi warps but no nearby, massive companions (though, in the cases 
of NGC3432 and NGC3079, some associated or nearby low surface brightness or dwarf 
satellites). These disk systems resemble NGC5907, which at the time of this pilot survey 
had already shown to have a tidal loop by Shang et al. (1998), and, in the case of NGC 3044, 
a minor merger was already hypothesized (Lee & Irwin 1997). 

The observations were made with the Mosaic camera during several nights (UT 2000- 

03-31 to 2000-04-04) using the BATC 9 filter — essentially a narrow version of the Cousins 
R band — used by Shang et al. (1998) in their imaging of NGC 5907; the filter was kindly 
loaned by Rogier Windhorst. The Mosaic camera is made of eight 2k x 2k individual chips 
and provides a field of view (FOV) of 59' at a pixel scale of 0.43"pixel~^. The conditions 
for this KPNO run were photometric, with a typical seeing of ~ 1.3 ". Using standard 
IRAF procedures for the over-scan/bias correction and flat- fielding, we combined the final 
image (see upper panel of Fig.l) of NGC 4013 from the sum of four 1200 sec exposures, eight 
1800 sec exposures, and one 900 sec exposure. 

A very faint arc is clearly identified in this image, situated at ~ 6' northeast of the 
center of the galaxy in the region of the extreme of northern plume of the gaseous warp 
(Bottema 1996, see Sec. 4). However, as may be seen, the achievable fiat- fielding and surface 
brightness limit was not adequate to make definitive conclusions about the nature of this 
structure. 



2.2. Isaac Newton Telescope 

To confirm and better trace the extent of this structure, follow-up observations were 
obtained with the 2.5m Isaac Newton Telescope (INT) at La Palma using the Wide Field 
Camera (WFC) at f/3.29. This instrument holds four 4096' x 2048' pixels EEV CCDs with 
a pixel size of 0.332", providing a total field of about 35'x35'. Images were acquired during 
two service nights (UT 2003-03-27 and 2003-04-03) using the Sloan r' (#214) and Harris B 
(#191) filters. From the first night we have two 600 s r' band exposures and from the second 
night three additional 600 s r' band images plus four 1350 s B band exposures. 

The central CCD (chip 4) was large enough (FOV 22.8'xll.4') to include NGC 4013 and 
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its warp, so we reduced only chip 4 in isolation, using standard IRAF procedures for the over- 
scan/bias correction, flat-fielding, and combining of the individual images. The background 
on the final B band image turned out to be fiat and we obtained a sky value of 22.22 B- 
inag/n" from ellipse fits, using the fixed outer disk geometry, beyond the galaxy as described 
in Pohlen & Trujillo (2006). The resulting 3-sigma uncertainty on the sky value is below 
0.04%. However, the individual R band exposures suffered from scattered light causing large- 
scale variations in the background of the order of 2%. After careful masking of all sources, 
which were replaced by the mean of a linear interpolation along lines and columns, we used 
the IRAF imsurfit routine on a median filtered (~ 30" window) version of the masked and 
interpolated image to determine the background on each individual one. After subtracting 
this second order fit from the images they were combined to the final r' band version reducing 
the large-scale variations in the background to 0.2%. The two images used in this study have 
an equivalent exposure time of 50 min for the r' band and 90 min for the B band. Due to 
the unstable weather conditions (high humidity, sometimes cirrus) during the two nights we 
did not use the observed standard stars, but obtained our photometric calibration from the 
r' and g' band Sloan Digital Sky Survey (SDSS) (York et al. 2000) by means of aperture 
photometry. After geometrically mapping the two INT and two SDSS images to the same 
scale and center and applying an identically-sized, but conservative (i.e. large) mask on the 
disturbing, bright star close to the center of NGC4013, we obtained fiuxes in six concentric 
apertures (40"-240"diameter) on all four images. Following Smith et al. (2002) we converted 
the magnitudes measured on the SDSS images to Cousins R and Johnson B, and obtained 
the photometric zero-points for our INT images using a linear fit without a color term. The 
uncertainty is of the order of 0.01 mag. We get similar results by using aperture photometry 
of six isolated, bright, but unsaturated stars. All the following magnitudes are given in R 
and B, correcting for Galactic extinction, {Ab — 0.072, Ar — 0.044; Schlegel et al. 1998), 
but not for inclination. 

The final image of NGC 4013 is shown in the lower panel of Figure 1. Although the 
background is still an issue we can now more confidently trace and measure the low surface 
brightness arc. It is important to remark that all quoted surface brightness measurements 
in this paper (see Sec. 3) come from the data set obtained at this telescope, the only one 
where we have calibrated, deep enough photometric data. 

2.3. Black Bird Remote Observatory Telescope 

The evidence of a possible tidal stream around NGC 4013 from our previous data ob- 
tained with the KPNO telescope and INT motivates us to select this galaxy as a priority 
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target in our pilot survey of stellar tidal streams in nearby spiral galaxies using small tele- 
scopes, of which searching strategy was successful proven with the detection of very faint 
parts of the NGC 5907 tidal stream (the commissioning target of this survey; see Martmez- 
Delgado et al. 2008). 

With this purpose, we obtained very deep images of NGC4013 with the f/8.3 Ritchey- 
Chretien .5-meter telescope of the Black Bird Remote Observatory (BBRO) situated in the 
Sacramento Mountains (New Mexico, USA). The data was obtained during several dark 
sky observing runs between UT 2006-11-06 and UT 2006-12-28. We used a Santa Barbara 
Instrument Group (SBIG) STL- 11000 CCD camera, which yields a large field of view (27.7' 
X 18.2') and a plate scale of 0.45"pixel^^ The data consists of multiple deep exposures 
through non-infrared clear luminance (CL; 3500< A < 8500), red, green and blue filters from 
the SBIG Custom Scientific filters set. Table 1 gives a summary of the images collected for 
this project. Column 5 refers to the total exposure time of the co-added images for each 
filter obtained with each run. This yields a total exposure time per filter of 5400s (red), 
3240s (green), 6480s (blue) and 39600s (CL). 

A master dark and bias frame was created by combining 10 dark sub-exposures each 
produced at the same exposure length and camera temperature settings used for the lumi- 
nance and the filtered exposures. A master flat was produced by combining 10 separate sky 
flat exposures for each fllter. The science data sub-exposures were reduced using standard 
procedures for bias correction and flat flelding. The data captured through the CL fllter was 
median combined to produce a master CL-fllter data set. The red, green and blue flltered 
exposures were separately combined (using the median technique) to produce individual red, 
green and blue master image that represented the total exposure time through each color 
filter. The resulting data sets for each color filter were subsequently summed to produce a 
"synthetic" luminance image that represented the total exposure time accumulated through 
all of the color filters, that was subsequently summed combined with the CL-filter image . 
The resulting final image thus represented the sum of all available CCD exposures collected 
in this project, with an accumulated exposure time of 15.2 hours (including 11 hours through 
a clear luminance filter). 

The dynamic range of the final image was very large, because of the presence of a very 
bright, galactic disk surrounded by very diffuse tidal stream. To optimize the contrast of the 
faint structures detected around NGC 4013, we used an histogram equalization of the image 
by means of a non-linear transfer function^. This well-known image processing tecnique 
employs a monotonic, non-linear mapping which re-assigns the intensity values of pixels in 



^see http://http://homepages. inf.eci.ac.uk/rbf/HIPR2/histeq.htm and references within. 
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the input image such that the output image contains an uniform distribution of intensities. 
This yields an effective method to suppress the bright regions of the image, intensifying the 
fainter parts of the stream. In addition, to increase the signal-to-noise of the faint structures 
surrounding NGC4013, image noise effects were reduced by the application of a Gaussian 
filter (Davies 1990). An illustrative example of the effect of these processes is given in Fig. 2. 

The final image of NGC 4013 was obtained using an iterative process that involves 
several passes of a S-shaped histogram operator to the total image and the application of a 
three-pixel diameter low-pass Gaussian filter. This approach optimized contrast and detail 
throughout the tonal range of the image and permitted the faint structure surrounding 
NGC 4013 to be rendered distinguishable in the final image. The resulting image is shown 
in Figure 3 . We have added the labels A through G to identify some photometric feature 
we discuss in the following Section. 

As it is discussed in Martinez-Delgado et al. (2008), the search strategy of this survey 
was designed to obtain the position and morphology of very faint structures around the spiral 
galaxy, with the aim of providing a guide for pencil-beam, follow-up photometric surveys in 
state-of-art telescopes. Because of this approach is mainly based on the use of a very broad- 
band, uncalibrated luminance filter, the BBRO telescope data is only used in this paper to 
discuss the overall morphology of the possible stream. From our previous experience (e.g. 
NGC 5907: Martinez-Delgado ct al. 2008; Messier 94: Trujillo et al. 2008, in preparation), 
we estimate that the final image should include diffuse light structures as faint as /xr ~ 29 
mag/ . 

3. THE TIDAL STREAM(S) OF NGC 4013 

In addition to the well known peanut-shaped bulge (van der Kruit & Searle 1982, feature 
A in our Fig. 3) with a full vertical extent of about 2kpc (30"), a contour plot of the INT 
image (Fig. 4) reveals a huge, at least 13kpc (> 3') sized box-shaped outer "stellar halo" 
(feature B in Fig. 3 and also obvious in the images in Fig.l). Furthermore, to the East, a 
striking horseshoe-like structure (feature C in Fig. 5) is visible in the BBRO image, starting 
at the northeast corner of the box-shaped halo and reaching out to at least 26kpc (6') from 
the center. This feature, that corresponds to the arc detected in the KPNO and INT images 
(see Sec. 2 and Fig.l), is now evident as loop-like. The obvious hole in the light distribution 
clearly shows that this feature is not simply a stellar warp, but an actual "ring-likestructure 
seen obliquely, and of the type produced by tidal streams (see Sec. 4). The loop apparently 
enters the galaxy again at the southeast corner of the box (labeled F in Fig. 3). 
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The BBRO and INT images also show a pair of wings jutting to the west and southwest 
(features D and E) and straddhng this side of the NGC 4013 disk, squaring off the most 
diffuse galaxy light distribution at large radii. In the diffuse light visible in the northwest 
edge of the galaxy (feature E), there is some evidence for a shorter, coherent loop feature 
visible in a high-contrasted version of the BBRO image (Fig. 5). This is consistent with 
the presence of a second arm of debris (as expected in theoretical simulations, see Sec. 4), 
which is probably related to the northeast loop (feature C). However, deeper images are 
necessary to confirm this hypothesis. The more plume-like feature visible on the southwest 
side (feature D in Fig. 3) extends ~ 17 kpc (4') out. Unfortunately, the clarity of this feature 
(also barely visible in the INT image in Fig. 1) is somewhat hindered by the presence of two 
bright stars. 

Together, these wings, and the obvious loop to the northeast could be seen as parts 
of one continuous ribbon-bow or pretzel-like structure surrounding the disk but with arms 
that cross each other at the galaxy center (at least in projection). Evidently, these particular 
features of the extended light distribution are the result of a single, coherent tidal disruption 
event, which is creating a low- inclination stellar stream around the disk of NGC 4013. An 
example of how such a tidal structure might be formed is shown in a tidal disruption model 
(see Figure 6 and discussion in Section 4). It is worth noting that we do not observe any 
bright spot that could be identified with a remaining, intact dwarf galaxy core among any 
of the low-surface brightness features. Therefore, our images do not provide any evidence 
on the final fate of the progenitor galaxy, which could be hidden inside the disk or could be 
completely disrupted by now. 

The surface brightness of the described structures is very close to the remaining in- 
homogeneities in the background of our CCD INT images, specially on the r'-band image, 
which makes it extremely difficult to determine accurate surface brightnesses. Using a set of 
small ~ 4"aA" boxes, strategically placed along the stream by avoiding stars and any residual 
obvious brightness condensation, we estimated the mean surface brightness of the stream 
on the east side (feature C) to be about /ib = 28.61q'4 /^R — 27.0lo'2. The large errors 
arc due to the uncertainty of ~ 0.1% in the local sky determination (a set of similar small 
boxes placed outside, but close to the stream) and are of the same order as the brightness 
variation along the stream. This yields a (B — R) color of 1.61q;4, which is a typically red 
value found for SO galaxies (e.g. Barway et al. 2005), and would be consistent with the redder 
dwarfs of the Local Group (Mateo 1998). Although the color of the loop is about 0.6 mag 
redder than the outer (radially/ vertically) parts of the disk/halo of NGC 4013 and only close 
to the mid-plane of the galaxy (where the dust lane certainly reddens the intrinsic colors), 
we find similar colors as for the stream within the attached large error. Therefore, deeper 
broad-band observations are needed to undoubtedly conclude whether its stellar population 
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is not composed of the same mix of stars found in the outer region of NGC4013's disk, an 
additional evidence that this horseshoe-hke feature did not originate in the disk itself. 

The width of this giant stream is variable along its path, but this may be due to a 
projection effect. A rough estimate for the full-width-half-maximum (FWHM) of the stream 
1.0 ±0.3 kpc (~ 14") is taken at the most distant part on the northeast side of feature C in 
Fig. 3. 

Figure 7 shows the position of the possible tidal stream with respect to the prominent 
Hi gas warp (Bottema et al. 1987; Bottema 1996). This comparison clearly shows that there 
is no Hi gas associated with the detected giant stellar feature (such as is the case in NGC 
3310; Wehner & Gallagher 2005) on the northeast side (feature C). On this side, the projected 
path of the stream shows a very similar inclination to the warped gas disk, but is clearly 
outside and almost enclosing (in projection) the measured Hi gas 

Interestingly, our high-contrasted version of the BBRO image (Fig. 5) reveals the trace 
of a spike- like feature in the northeast edge of the galaxy disk (feature G), with a position and 
orientation with respect to the galactic plane that is consistent with the stellar counterpart 
of the prodigious gaseous warp. It is clear, however, that (at our surface brightness limit) 
this stellar component warp is significant less extended than the Hi -component displayed 
in Fig. 7. On the southeast side the extended optical feature D (see Fig. 3) seems to be 
associated to the southern H I gas plume of the warp, although on this side we are severely 
hampered by the two bright foreground stars. However, the clear separation of the H I gas 
and the optical light on the other side of the galaxy let us reject the direct association of 
these structures as an extended, low luminosity stellar counterpart of the warped H I disk of 
NGC 4013. 



4. DISCUSSION 

Our deep images of NGC 4013 show a complex stellar halo, including a coherent, stream- 
like structure inclined ~ 25 deg with respect to the galaxy disk. We have independently 
confirmed the presence of these low surface brightness features with three data sets from three 
different telescopes. In addition, if one knows what to look for, the loop is also barely visible 
on the Digitized Sky Survey 2 (DSS2) blue plate and on the archival 4.5/i Spitzcr Infrared 
Telescope Facility image. The existence of the loop structure is therefore without question. 
The discovered structures are evidently debris material produced by the tidal disruption of a 
low-mass galaxy companion to NGC 4013: Such features are expected from tidal disruption of 
a companion galaxy (Johnston et al. 2001) and similar features are also observed in NGC 5907 
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(Shang et al. 1998; Martmez-Delgado et al. 2008), NGC3310 (Wehner & Gallagher 2005), 
and a handful of other external galaxies (summarized in Martmez-Delgado et al. 2008, Table 

1)- 

The sky-projected geometry of this structure is in fact similar to the edge-on, external 
perspective of a tidal stream that results from the tidal disruption of a single dwarf galaxy 
satellite moving on a low-inclined orbit predicted by different theoretical simulations (e.g. 
the Monoceros tidal stream: Penarrubia et al. 2005;Canis Major over-density: Martin et 
al. 2005, their Fig. 14). For comparison purpose, we will use here as template the N-body 
model constructed by Penarrubia et al. (2005) to reproduce the geometry and kinematics 
of the Monoceros tidal stream. This model simulates the accretion of a low-mass satellite 
(6 xlO^ Mq) on a nearly circular (e=0.10±0.05), low inchned (i ~ 25°) orbit in the last 
3 Gyr and provides a remarkable good match to the projected, edge-on view of the stream 
wrapping NGC 4013 (see Fig. 6). Considering that NGC 4013 has a smaller disk scale-length, 
Hr = 2.8 kpc, and a lower maximum circular velocity, K,max = 195 kms~^ (Bottema 1996) 
than the Milky Way, we can apply dynamical considerations to estimate the stream age 
from the N-body model, so that Age/AgeMon = hn/hR^MW x K^max/K.max.Afw- For the 
Milky Way we have that Hr = 3.5 kpc and Vc^max = 220 km/s, thus Age/Age^^^„ ^ 0.9. 
Since some pieces of the Monoceros stream are Age^^^^ = 3 Gyr old, we estimate that the 
stream surrounding NGC 4013 may be as old as ~ 2.8 Gyr. However, although the model 
provides a reasonable qualitative resemblance to the morphology of NGC 4013, it cannot 
fully explain the observed structure. For example, while the observed tidal structure in Fig. 
3 is clearly more prominent on its Eastern side of the galaxy, the model shows loops of 
roughly equal brightness on both sides (as well as various other features that not detected in 
our deep image). Unfortunately, with no kinematic data and only the projected geometry of 
the stream as constraint, the orbit and mass of the progenitor system may be degenerated, 
and in that respect the estimates presented here must be considered illustrative only. 

Although the progenitor system of the Monoceros and NGC 4013 tidal streams might 
have shared similar orbits and initial masses, our estimates reveal that the surface brightness 
of the latter is approximately ~2-3 magnitudes higher than the value reported for this low 
galactic latitude tidal debris candidate (/^b ~ 30.0mag/n": J. de Jong, private communica- 
tion). A possible explanation might be that the progenitor of the NGC 4013 stream was more 
luminous than that of the Monoceros stream. Since both models have the same initial mass, 
this hypothesis would imply that the mass-to-hght ratio of the NGC 4013 stream progenitor 
was considerably lower than that of the Monoceros stream. However, the N-body models 
used by Peiiarrubia et al. (2005) assume that the stellar and dark matter components of 
the progenitor system have the same spatial distribution (i.e. mass- follows- light models). In 
a ACDM cosmogony, stars only populate the inner-most regions of dwarf galaxies (Strigari 
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et al. 2007, Penarrubia et al. 2008a) and can only be tidally stripped after most of the 
dark matter halo beyond the luminous radius has been lost (Penarrubia et al. 2008b), which 
introduces an additional parameter: the spatial segregation of the stellar component with 
respect to the surrounding dark matter halo of the progenitor galaxy, which is fundamental 
to determine the survival time of dwarf galaxies and the properties of their associated tidal 
streams. For example, the more deeply embedded within the halo a stellar component is, 
the thinner, brighter and colder will be its associated tidal stream. Therefore, for ACDM- 
motivated models, the halo size, the stellar spatial segregation and the total luminosity of a 
dwarf galaxy are all parameters that influence the resulting surface brightness of an associ- 
ated tidal stream but which, unfortunately, cannot be constrained in absence of additional 
data (e.g. kinematics). 

Whereas the detection of the stream around NGC4013 cannot shed hght on the nature 
of the Milky Way stellar over-densities (including on the controversial origin of the Monoceros 
stream; see Sec. 1), it does suggest that the accretion of satellite galaxies on low-inclination, 
low-eccentric orbit may be relatively common during the hierarchical formation of spiral 
galaxies. In fact, their occurrence in the known sample of stellar tidal streams in the Local 
Volume (D<15Mpc; see Martmez-Delgado et al. 2008; Pohlen et al. 2004) is >20%. In this 
context, the presence of giant disks around spiral galaxies that extend out to several scale- 
lengths of the inner disk (e.g. the extended disk found by Ibata et al. (2005) around M31 is 
detectable out to a projected radius of 80 kpc with the present instrumentation) might also 
be signature of the accretion of massive satellite galaxies moving on low-inclination, low- 
eccentricity orbits (Penarrubia et al. 2006), adding support to the scenario that the disks of 
spiral galaxies may grow inside out as a result of this type of accretion events. 

The detection of galactic warps in spiral galaxies that present signatures of recent accre- 
tion events is quickly increasing with the availability of wide-field, deep photometric surveys. 
Although the origin of galactic warps is still unknown, it was early proposed that the gravi- 
tational perturbations induced by satellite galaxies may cause the formation of warps (Burke 
1986, although see Hunter & Toomrc 1969 for counter- arguments). Since those early studies, 
several alternative scenarios have been proposed in the literature. For example, it has been 
shown that the interaction of a stellar disk with the surrounding dark matter halo may also 
induces the formation of warps (Sparke & Casertano 1988; Binney et al. 1998), as well as 
bending instabilities (Revaz & Pfenniger 2004), intergalactic accretion flows onto the disk 
(Lopez-Corredoira et al. 2002) and cosmic in-fall (Shen & Sellwood 2006). 

A crucial fact that has been often invoked to dismiss the tidal origin of galactic warps was 
the existence of warped spiral galaxies in apparent isolation (Sancisi 1976). NGC4013, with 
its prominent, and rather symmetrical H I warp, along with a really unassailable indication of 
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a past merger event, provides an important counter-example to the above argument; indeed, 
NGC4013 may well be the Rosetta Stone for warp theories, since morphologically the warp 
and the merger debris seem so closely aligned. 

With NGC 5907 and NGC 4013 we have now two examples of apparently isolated galax- 
ies with significantly warped gaseous disks, which at the same time show evidence for an 
ongoing tidal disruption of a dwarf companion. The connection between the warp of the 
MW and its satellite galaxies has been already explored in some recent studies. For exam- 
ple, Bailin (2003) and Weinberg & Blitz (2006) study the effects of the Sagittarius dwarf and 
the Magellanic Clouds on the MW disk, respectively. Interestingly, the existence of warped 
galaxies that have suffered a recent accretion of satellite galaxies on low-inclination, low- 
eccentricity orbits (like the MW and NGC 4013) suggests that these kind of accretion events 
may induce the formation of warps more efficiently than the perturbations from satellite 
galaxies moving on highly eccentric, nearly-polar orbits. This is a possibility worth to be 
further investigated. 

With the growing numbers of examples showing a connection between warped disks and 
evidence of mergers, even for disk systems with no obvious companions, we may be close to 
resolving the mystery of warps, and in a way be consistent with prevailing ACDM models 
of galaxy formation. 
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Tablc 1. JOURNAL OF OBSERVATIONS (BLACK BIRD OBSERVATORY) 
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Fig. 1. — NGC4013: Smoothed and enhanced versions of the KPNO 0.9-m image (upper 
panel) and the INT B band image (lower panel) highhghting the low-surface brightness 
features. The arrows of lengths 2' in the lower left corner give the size and orientation. Both 
images are displayed in logarithm scale. 
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Fig. 2. — An illustrative example of the effect of histogram equalization and Gaussian filtering 
techniques in the resulting final image obtained from the sum of all the exposures taken with 
the BBRO: A) the image with linear-stretching (and thus non-proccesed) ; B) the image 
after applying an histogram equalization with a non-linear function; and C) the effect of a 
Gaussian filter of 5 pixels of radius in the summed image displayed in panel A. The blurrer 
aspect of this image with respect to the final version displayed in Fig. 3 is due the larger 
filter widthness used here to enhancent the effect of the filter with illustrative purposes. The 
field of view of these images is 28' x 18'. 
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Fig. 3. — Image of NGC4013 obtained with the BBRO 20-inch telescope. The total exposure 
time of the original image was 13.7 hours (including 11 hours in clear-luminance filter) and 
was noise-filtered by applying a Gaussian filter. The image has dimensions of ~ 16.8' x 12.8', 
which, at the distance of NGC 4013 is ~ 71 x 54 kpc. Identified photometric features labeled 
A-G are discussed in Sec. 3. For illustrative purpose, a color image of the NGC 4013 obtained 
with the same telescope has been superimposed on the saturated disk region of the galaxy. 
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Fig. 4. — Isaac Newton Telescope B band contour map of NGC4013 highlighting the box- 
shaped outer halo. Contours are at levels 17.0-27.0 B-mag arcsee"^ in equidistant steps of 
0.5 mag. Towards the outer parts we smoothed the contours using three increasing levels of 
median filtering (indicated by the three different colors). 
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Fig. 5. — A high- stretched version of the luminance clear-filter image of NGC4013 obtained 
with the BBRO telescope. This shows evidence of a possible second arm of debris within 
the diffuse light wind to the northwest edge of the NGC4013 disk (feature E in Fig. 3). In 
addition, the spike-like feature visible in the northeast edge of the disk could be related to 
the stellar component of the prodigious gaseous warp of the galaxy. North is up and East 
to the left. 
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Fig. 7. — Overlay of the Hi contours from Bottema (1996) (his Fig.l) on our BBO image of 
NGC4013. The lowest Hi contour level shown is about 1x10^° H-atoms cm~^. The beam is 
given in the lower right corner of the Hi insert. The scale of this image is ~4.2 kpc/arcmin. 



